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We study the decays B + — > n c K + and B° — > n c K°, where the n c is reconstructed in the K®K ± ir^ : 
and K + K~tt° decay modes. Results are based on a sample of 86 million BB pairs collected with 
the BABAR detector at the SLAC e + e~ B Factory. We measure the branching fractions B{B + — ► 
ri c K + ) = (1.34 ±0.09 ±0.13 ±0.41) x 10~ 3 and B(B° -> ri c K°) = (1.18 ±0.16 ±0.13 ±0.37) x 10~ 3 , 
where the first error is statistical, the second is systematic, and the third reflects the r\ c branching 
fraction uncertainty. In addition, we search for B — > r\ c K events with r\ c — > 2(K + K~) and n c — * 4>4> 
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and determine the ij c decay branching fraction ratios B(r\ c — ► 2{K + K ))/B(rj c —* KKtv) — (2.3 ± 
0.7 ±0.6) x 10~ 2 and B{ij c -Kj><t>)/B(r) c -> KKtv) = (5.5 ± 1.4 ± 0.5) x 10 -2 . 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 

The decay B — > j] e K is used to measure sin2/3 Q, @ , 
but is interesting dynamically as well. The ratio of its 
decay rate to that of B — > J/ipK reflects the underly- 
ing strong dynamics and can be used to check models of 
heavy quark systems [E El IE IE • The strong decay 
should be isospin invariant, an expectation that can be 
checked and then used to combine results for higher preci- 
sion. It is therefore interesting to measure accurately the 
branching fractions for B° — > rj c K and B + — > r/ c K + 0. 

We use data collected with the BABAR detector at the 
PEP-II energy-asymmetric e + e _ storage rings. The data 
sample contains 86.1 x 10 6 BB pairs, corresponding to 
an integrated luminosity of 79.4 fb _1 taken at a center- 
of-mass energy equivalent to the mass of the T(4B) reso- 
nance. An additional 9.6 fb _1 of data, collected 40 MeV 
below the resonance, is used to study the background 
from light quark and cc production. 

A detailed description of the BABAR detector can be 
found elsewhere 9]; only detector components relevant 
to this analysis are mentioned here. Charged-particle 
trajectories are measured by a five-layer double-sided sil- 
icon vertex tracker (SVT) and a 40-layer drift cham- 
ber (DCH), operating in the field of a 1.5-T solenoid. 
Charged particles are identified by combining measure- 
ments of ionization energy loss (dE/dx) in the DCH and 
SVT with angular information from a detector of in- 
ternally reflected Cherenkov light (DIRC). Photons are 
identified as isolated electromagnetic showers in a CsI(Tl) 
electromagnetic calorimeter. 

In this analysis, r\ c mesons are reconstructed in the 
K^K^, K+K-tt , 2{K+K-) and (jxf) decay modes. 
Candidates for K® are identified through the decay — > 
7r + 7r~, (j> candidates through <f> — ► K + K~ and tt° can- 
didates through 7T° — > 77. Note that ?y c decays to 
2(K + K~) include both non-resonant and resonant {^>^>, 
<j)K + K~) components, so we expect a partial overlap of 
the 2{K + K~) and cfxf> samples. 

We require that charged tracks, other than those used 
to reconstruct if" — * 7r + 7r~ candidates, have a minimum 
transverse momentum of 0.1 GeV/c, and that they origi- 
nate from the interaction point to within 10 cm along the 
beam direction and 1.5 cm in the transverse plane. The 
"fast" kaon candidate at the two-body B + decay vertex 
is required to have at least 12 hits in the drift chamber 
and to have a momentum in the T(AS) rest frame larger 
than 1.5 GeV/c. 
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The cuts used to select K + , K®, <p and ir° candidates 
from rj c decays are described below. They are optimized 
to maximize the statistical sensitivity of the signal, de- 
fined as S/^JS + B, with S and B being the estimated 
numbers of signal and combinatorial-background events. 

All charged-kaon candidates are required to have mo- 
mentum greater than 250 MeV/c and a polar angle be- 
tween 0.35 and 2.54 rad with respect to the detector 
axis, to restrict them to a fiducial region where the par- 
ticle identification performance can be determined with 
small uncertainty. Kaon identification is based on a neu- 
ral network (NN) algorithm that combines information 
from the DCH, the SVT, and the DIRC. Particle iden- 
tification criteria are crucial for background suppression, 
especially for r\ c — ► 2(K + K~) decays. In this channel, 
three of the four kaons must pass a tight cut on the NN 
output variable. Less restrictive requirements on the NN 
signature are used for identifying the fourth kaon from 
■q c — > 2(K + K~) candidates, the charged kaons in the 
other rj c decay modes, and the fast kaon from B + de- 
cays. The kaon-identification efficiency depends on the 
momentum and polar angle of the track, as well as on 
the chosen NN cut. For the tightest selection above, the 
average kaon efficiency exceeds 85%; the corresponding 
pion-rejection efficiency is about 98%. 

We reconstruct K® candidates from pairs of oppositely 
charged tracks fitted to a common vertex. We require 
the candidate from the B (rj c ) decay to have a recon- 
structed mass within 13 (16) MeV/c 2 of the K® mass [Toj . 
Furthermore, the cosine of the opening angle between the 
flight direction and the momentum vector of the can- 
didate is required to be larger than 0.9995 (0.9930), and 
the flight distance from the B vertex larger than four 
times its error. 

We reconstruct (f> candidates from pairs of oppositely 
charged kaons with an invariant mass within 14 MeV/c 2 
of the <f> mass [Toj . 

We use pairs of photons to reconstruct ir° — > 77 can- 
didates, requiring a minimum energy of 120 MeV for one 
photon and 80 MeV for the other. The reconstructed 
77 mass is required to lie within 18 MeV/c 2 of the 7T° 
mass 10]. 

We reconstruct rj c candidates by fitting the appropriate 
combination of charged tracks, K®, cf>, or ir° candidates 
to a common vertex. Neutral or charged B candidates 
are formed from reconstructed ?y c and K® or K + candi- 
dates. In reconstructing the B decay chain, the measured 
momentum vector of each intermediate particle is deter- 
mined by refitting the momenta of its daughters, con- 
straining the mass to the nominal mass of the particle 
and requiring that the decay products originate from a 
common point. In the case of the rj c , only the geometrical 
vertex constraint is applied because of the large intrinsic 
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width of the resonance. Charmonium candidates are ac- 
cepted if they have an invariant mass between 2.7 and 3.3 
GeV/c 2 . Note that this procedure also reconstructs J/tp 
decays, which are used to measure the mass resolution 
and for other cross-checks. 

We use a Fisher discriminant to suppress e + e~ — > qq 
background processes. Our Fisher discriminant is a lin- 
ear combination of 18 variables, the most important of 
which are the normalized second Fox- Wolfram moment 
and the angle between the thrust axis of the B candidate 
and that of the rest of the event. Also contributing are 
the energy flow in nine 10° polar angle intervals coaxial 
around the r\ c direction in the center-of-mass frame |l l| , 
the polar angles of the B candidate and of the overall 
thrust axis, and other event-shape variables that distin- 
guish between BB events and continuum background. 
The discriminant is tuned on simulated signal events and 
on off-resonance data to achieve maximum separation 
between signal and continuum background. Fisher co- 
efficients arc determined individually for each rj c decay 
mode, and threshold values are set as part of the cut- 
optimization procedure described earlier. 

We select B candidates using two nearly independent 
kinematic variables: toes, the beam-energy-substituted 
mass, and AE, the difference between the energy of the 
B and the beam energy in the center-of-mass frame 0. 
The Toes resolution is 2.6 MeV/c 2 , dominated by the 
beam-energy spread. The AE resolution varies from 
15 to 28 MeV, depending on the rj c decay mode. Sig- 
nal events are expected to have ities close to the B 
mass and AE close to zero. Our selection requires 
5.2 < toes < 5.3 GeV/c 2 . After all the cuts listed above, 
10-25% of the selected events, depending on the rj c decay 
channel, contain more than one B — > rj c K candidate in 
a AE window ±250 MeV wide; we then retain only the 
candidate with the smallest value of |Ai?|. We have veri- 
fied with simulated events that this procedure selects the 
correct candidate in 90-98% of the cases, and that it does 
not bias the measurement. Finally, we require candidates 
to lie within an optimized interval of AE that varies from 
±35 to ±70 MeV, depending on the decay mode. 

Events surviving the full selection chain originate from 
four different sources: B — ► r\ c K decays, i.e., the signal; 
B — > J/ipK decays, with the J/ip decaying into the same 
final state as the r\ c \ a combinatorial background, aris- 
ing from random track combinations in continuum and 
in BB final states; and a background component from 
other B decays to the same final state particles as the 
B — ► rj c K decay mode under consideration. The last 
background component can contribute events that clus- 
ter at the signal peak in Toes and AE and is therefore 
termed "peaking background" . 

Examples of peaking background for B — > rj c K (t] c — > 
K^K^) are B+ -> K*~ K+ K+ (K*~ -> K° s n ) or 
B° -> K*°K°K° (K*° -► K+tt ). In the particu- 
lar case of the decay B + — > rj c K + (r\ c — > K®K x^), 
another important source of peaking background comes 
from B+ -> D°K° S K + (D° -> K+tt'). For this B de- 



cay mode therefore, candidates with a K + ix~ invariant 
mass within 15 MeV (3a) of the D° mass are explicitly 
vetoed. Other processes, such as nonresonant B decays 
to the selected final state, whose branching fractions are 
not well-known, can also contribute. The mass mx of the 
system recoiling against the fast kaon is used to separate 
B — > rjc K and B — > J/ipK events, which peak at the 
mass of the corresponding charmonium system, from the 
peaking (in toes and AE) background, which is expected 
to exhibit a linear dependence on mx- This assumption 
is verified with large samples of simulated BB events. 
These studies also show that inclusive B decays into r\ c 
and potential cross-feed among different r\ c decay modes 
are negligible after the event selection. 

The number of signal events is determined from an un- 
binned maximum-likelihood fit to the joint Toes and mx 
distribution. Four hypotheses are considered to build the 
2-D likelihood function: B — > r\ c K signal, modeled by the 
product of a Gaussian resolution function in toes and of 
a non-relativistic Breit-Wigner function convoluted with 
a Gaussian resolution function in mx ; B — ► J/tp K com- 
ponent, given by the product of Gaussian resolution func- 
tions in toes and mx] combinatorial background, mod- 
eled by an "ARGUS" endpoint function in toes 03 > and 
a linear function in mx] and peaking background, de- 
scribed by a function linear in mx and Gaussian in Toes- 
The widths of the mx and Toes resolution functions, and 
the mean value of the Toes distribution are free parame- 
ters common to the r\ c and J/ip probability density func- 
tions (p.d.f.). The latter two parameters also determine 
the toes dependence of the peaking-background p.d.f., re- 
flecting the evidence that this background is dominated 
by B decays to the same final states as the signal. We set 
the r] c and J/ip masses to their world-average values ^3 > 
the endpoint of the combinatorial background function 
to 5.29 GeV/c 2 , and the rj c width to the value recently 
measured by BABAR ^|. All other parameters and the 
number of events in the different components are deter- 
mined by the fit, which is performed separately for each 
decay channel. 

For B — ► rj e K modes with r\ c — > KKir, candidates 
are weighted to take into account small efficiency varia- 
tions across the rj c Dalitz plot. The weighting procedure 
compensates for any resonant structure in r/ c three-body 
decays unaccounted for by the simulated phase-space dis- 
tribution, which is uniform over the Dalitz plot. Since all 
weights are close to one, they do not affect the shape of 
the different components and have only a marginal in- 
fluence (0.6-4%) on the fitted event yield. Samples of 
simulated events are used to verify that the likelihood fit 
is unbiased. 

The measured i] c signal yields are reported in Table [I] 
We observe a significant signal in all modes with the 
exception of B° with rj c decaying into 2(K + K~) and 
(pip. The toes and mx distributions of B + — > ij c K + 
candidates are shown in Fig. ^ In the largest samples 
(r) c — > K®K ir^) we can determine the r/ c width T(rj c ) 
from a simultaneous fit to neutral and charged B data. 
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TABLE I: Number of B — > rj c K events and statistical signif- 
icance <S, denned as ^/21og(£ max /£o) where £max/£-o is the 
likelihood ratio of the fitted maximum over the null hypothe- 
sis. The first error is statistical, the second is the systematic 
uncertainty associated with the fitting procedure. 



Mode 


?7c yield 


<S 


B+ Vc K + 






r) c -> KIR*** 


306.4 ± 24.4 ± 14.0 


20.5 


rj c -> K+K'n° 


136.8 ± 17.5 ± 9.3 


11.7 


Vc -» 2{K+K~) 


26.2 ± 8.4 ± 4.5 


4.5 


r/ c -» 00 


19.1 ± 4.9 ± 0.6 


6.6 










79.4 ± 12. 7± 4.3 


9.7 


40.9 ± 9.5 ± 2.7 


6.2 


r lc ->2(K+K-) 


3.9 ± 3.7 ± 1.6 


1.4 


r/ c — ► 00 


3.0 ± 1.7 ± 0.1 


3.6 



We find T(t] c ) = 39.7 ± 6.6 MeV/c 2 , where the error is 
statistical only, consistent with the BABAR measurement, 
r(?7 c ) = 34.3 ± 2.3 ± 0.9 MeV/c 2 0. 

The systematic uncertainty associated with the fitted 
signal yield includes three components: the uncertainty 
in the fixed parameters, the uncertainty associated with 
the Dalitz weighting procedure, and the uncertainty as- 
sociated with the p.d.f. models. The first component is 
evaluated by varying each fixed parameter, one at a time, 
by one standard deviation and repeating the fit. This 
component is dominated by the uncertainty on r(77 c ) (0- 
3% fractional uncertainty in B, depending on the mode). 
For the second component, the fit is repeated without 
applying the Dalitz-correction procedure; half of the dif- 
ference on the rj c signal yield (0-2% in B) is conserva- 
tively assigned as the corresponding systematic uncer- 
tainty. The last component is dominated by the uncer- 
tainty in the peaking background model. This error is 
evaluated by varying the assumed mx dependence from a 
first- to a second-order polynomial; it typically amounts 
to 4%, and exceeds 10% only for the r/ c -> 2(K + K~) 
modes. The error associated with the mx resolution 
function model (0-5%) is estimated by using, instead of 
a single Gaussian function fitted to the data, double- 
Gaussian resolution functions fitted to each simulated 
signal sample. 

Efficiencies are computed with simulated signal events 
that are reconstructed and selected using the same pro- 
cedure as for the data, including the yield-extraction 
fit. We apply small corrections, determined from data, 
to the efficiency calculation to account for the overesti- 
mation of the tracking and particle-identification perfor- 
mance, and of the 7r° and reconstruction efficiencies. 
A systematic uncertainty is assigned to each correction 
to account for the limited size and purity of the control 
sample used in computing that correction. For example, 
for the fast kaon identification, we correct the Simula- 
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FIG. 1: Distributions of tties (left) and mx (right) for charged 
B candidates. The tties distributions displayed here are re- 
stricted to the 2.90 < mx < 3.15 GeV/c 2 range; similarly, 
the mx distributions include only events in the toes sig- 
nal region (m ES > 5.27 GeV/c 2 ). Each pair corresponds 
to a different r\ c decay mode: (a, b) r/ c — * KsK ± tv^ 1 ; (c,d) 
77c K+K-ir ; (e,f) 77c -» 2(K + K~)- and (g,h) 77c 00. 
The fitted p.d.f. projections are shown as solid curves. In 
each plot, the dark grey region corresponds to the combina- 
torial background component, light grey highlights the peak- 
ing background, and the dotted line is the sum of the total 
background and of the J/ij) component. 



systematic uncertainty contributions associated with the 
sample size, the background subtraction, and the differ- 
ent kinematics of this decay chain compared to the two- 



7 



TABLE II: Efficiencies and reiative systematic uncertainties. 

B + r, c K + W~^J<1 

JjTg^Tr^ Jf + jy-7r 2{K+K~) 4><t> K%K^k^_ K + K~n° 2{K+K~) <j><p 

Signal efficiency 

0.213 0.124 0.155 0.194 0.184 0.126 0.147 0.170 



Source of uncertainty Relative uncertainty on signal efficiency(%) 



Monte Carlo statistics 


1.0 


1.5 


1.2 


1.1 


1.2 


1.3 


1.2 


1.1 


Tracking 


6.0 


3.4 


6.0 


6.0 


7.8 


5.2 


7.8 


7.8 


K® reconstruction 


3.0 








6.0 


3.0 


3.0 


3.0 


Particle identification 


3.9 


6.5 


12.1 


8.0 


1.5 


3.9 


9.4 


5.9 


7T° reconstruction 




5.0 








5.0 






Selection cuts 


2.7 


2.8 


1.7 


2.2 


3.3 


3.5 


3.1 


3.4 


Yield-extraction fit 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


Total uncertainty 


8.8 


9.8 


13.9 


10.7 


10.9 


9.9 


13.3 


11.2 



body B + — > rj c K + decay. Similarly, corrections affecting 
the 7T° reconstruction are calibrated using real and sim- 
ulated e + e~ — > t + t~ and multihadron samples. 

In addition, after all corrections, we compare our signal 
simulation to appropriate control samples with similar 
kinematics or final-state topology, in order to quantify 
the ability of the simulation to model the kinematic and 
event-shape variables used in the event selection. The 
small residual differences in the efficiencies at the cut 
value are assigned as systematic uncertainties affecting 
the selection procedure. 

Finally, we assign a systematic uncertainty to the yield- 
extraction fit by evaluating the influence of mixing back- 
ground events with simulated signal events. Values for 
the efficiencies, the corrections, and the corresponding 
systematic uncertainties are reported in Table ITT1 

The results on the products of the branching frac- 
tions for each mode are listed in Table ITTT1 We use the 



TABLE III: Measured branching-fraction products B(B — > 
r/ c K)xB(r] c — > X) (10 -6 ). The first error is statistical and 
the second is the total systematic uncertainty. 



Vc 


decay channel 


B + -> Vc K + 




Vc 


- K ii K-7T+ 


48.6 ± 3.9 ± 4.9 


42.6 ± 6.8 ± 5.2 


Vc 


-> K+K-it 


12.9 ± 1.7 ± 1.6 


11.1 ± 2.6 ± 1.3 


Vc 


-> 2{K + K~) 


2.0 ± 0.6 ± 0.4 


0.9 ± 0.9 ± 0.4 


Vc 




4.7 ± 1.2 ± 0.5 


2.4 ± 1.4 ± 0.3 



world-average values for the K® —> 7r + 7r~ , 7r° — > 77 and 
(f> — > K + K~ branching fractions ^3] and include their 
uncertainties in the systematic error. The systematic er- 
ror also comprises the uncertainties from the determina- 
tion of the number of BB pairs (1.1%), from the likeli- 
hood fit, and from the signal efficiency. We assume that 
the branching fraction of the T(4S) into BB is 100%, 
with an equal admixture of charged and neutral B final 
states. We do not include any additional uncertainty due 
to these assumptions. Possible interference effects be- 
tween the B — > rj c K signal and the peaking background 
are neglected. 



The decay amplitudes for i] c — > K + K~ir° and r\ c — > 
K°K~ir + are related by isospin symmetry. The ex- 
pected ratio of branching fractions, using the appropriate 
Clebsch-Gordon coefficients, is 0.25. Our measurements 
are consistent with this value for both the B + (0.27 ± 
0.04 ± 0.03) and the B° (0.26 ± 0.07 ± 0.03) sample. 
We therefore combine our results for these modes and 
use the world avera ge fo r the rj c — > KKtx branching frac- 
tion (0.055 ± 0.017113) to derive 

B{B + -> r lc K + ) = (1.34 ± 0.09 ± 0.13 ± 0.41) x 10~ 3 

B(B° 7] C K°) = (1.18 ± 0.16 ± 0.13 ± 0.37) x 10" 3 

where the first error is statistical, the second systematic, 
and the third due to the uncertainty on the r/ c — > KKt: 
branching fraction. In the combination we separate cor- 
related and uncorrelated uncertainties to weight the indi- 
vidual results and obtain the total systematic error. We 
also compute the ratio of neutral over charged B decays 
B(B° -> r] c K°)/B(B+ -> r) c K+) = 0.87 ± 0.13 ± 0.07, 
and, multiplying by the mean lifetime ratio t b + /t b o = 
1.085 ± 0.017 [l£|, we derive the ratio of partial widths 

T(B Q -> 7 lc K°)/T(B + -> T] C K+) = 0.94 ±0.14 ±0.08. 

To determine R K = T(B nJQ/T(B -> J/ipK), 
we use the BABAR measurements [lj| of the branching 
fractions, B(B+ -> J/^K+) = (10.1 ± 0.3 ± 0.5) x 10~ 4 
and B(B° -> J/^K°) = (8.5 ± 0.5 ± 0.6) x 10" 4 . Wc 
obtain 

R K (B + ) = 1.33 ± 0.10 ± 0.12 ± 0.41 

R K (B°) = 1.39 ± 0.20 ± 0.13 ± 0.43, 

where the first error is statistical, the second systematic, 
and the third due to rj c — > KKt: branching fraction. Our 
results agree with most predictions for Rk 1 which range 
from 0.9 to 2.3 0,0,|llEl3- 

The measured values of B(r] c — > 2{K + K~)) and 
BiVc —> 4"P) have higher uncertainties and therefore 
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these modes are not used for averages. We can ex- 
press our r] c — ► 2(K + K~) and r/ c — > 00 results in terms 
of ratios to the best-measured branching fractions of 
r; c — > KKtt, thereby cancelling all fully-correlated sys- 
tematic uncertainties. We average results on charged B 
decays and neutral B decays, taking into account corre- 
lations in the systematic uncertainties, to obtain B{rj c — > 
2(K+K-))/B(r) c -> KKtt) = (2.3±0.7±0.6) x 10" 2 and 
B(j] c -> H>)/B(ri c -> KKtt) = (5.5 ± 1.4 ± 0.5) x 10~ 2 . 
These results can be translated into rj c branching frac- 
tions: 

B{r) c -> 2(K + K~)) = (1.3 ± 0.4 ± 0.3 ± 0.4) x 10~ 3 

B(r) c -> 00) = (3.0 ± 0.8 ± 0.3 ± 0.9) x 10~ 3 , 

where the third error is due to the uncertainty of B(r\ c — > 
KKtt). Note that about half of the r/ c -> 2(K+K-) 
events are due to r\ c — > 00, — » K + K~ decays. Our 
measured branching fractions for r/ c — > 2(iT + ftr~) and 
Vc ~> <H> are consistent with recent results from Belle 
and BES ll5L llrl and are smaller than those of earlier 
experiments [lOj . 

As a cross-check, we can extract the branching fraction 
of J/tp decaying into the 2(K + K~) final state from the 
measured number of J/0 events in the appropriate J5 + 
and B° samples. Assuming the same efficiencies as for 
the B — > r] c K {rj c — + 2{K + K~)) processes and using 
the BABAR measurements of B(B + — ► J/t/jK + ) and of 
B(B° -> J/0 if ) 0, we obtain B(J/<tp -> 2(if+if-)) = 
(1.0 ± 0.5) x 10~ 3 (5+) and B(J/^ -> 2(iT+if-)) = 



(0.1 ± 0.2) x 10~ 3 (B°), where the error is statistical 
only. These results are consistent with the world average, 
B(J/ip -> 2(K+K-))= (0.7 ± 0.3 ) x 10" 3 [jjj. 

In summary, we have studied B — » ?7 c -ftr decays with 
?7 C decaying into i^i^TT^, K+K~tt°, 2[K + K~), and 
00. Using the first two decay channels, we have mea- 
sured the branching fractions B(B + — > t] c K + ) — (1.34 ± 
0.09 ± 0.13 ± 0.41) x 10~ 3 and B(B° -> r/ CJ ft: ) = (1.18 ± 
0.16 ± 0.13 ± 0.37) x 10~ 3 , which improve the statisti- 
cal precision of, and are in good agreement with, pre- 
vious measurements [B We have also measured 
the branching-fraction ratios B{r\ c — > 2(K + K~))/B(rj c — > 
KFtt) = (2.3±0.7±0.6)xlO- 2 andB(r; c -> cp(/>)/B(r) c -> 
-OTtt) = (5.5 ±1.4 ±0.5) x 10~ 2 , where ?? c -> 2{K+K-) 
includes ry c — * 00 events with — > K + K~ . The inferred 
branching fractions of ?7 C — > 2(X + i'C~) and ?y c — * 00 are 
in good agreement with recent results and smaller than 
suggested by earlier experiments. 
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